Plasmons in Sodium under Pressure: Increasing Departure from Nearly-Free-Electron Behavior 
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We have measured plasmon energies in Na under high pressure up to 43 GPa using inelastic x-ray scattering 
(IXS). The momentum-resolved results show clear deviations, growing with increasing pressure, from the pre- 
dictions for a nearly-free electron metal. Plasmon energy calculations based on first-principles electronic band 
structures and a quasi-classical plasmon model allow us to identify a pressure-induced increase in the electron- 
ion interaction and associated changes in the electronic band structure as the origin of these deviations, rather 
than effects of exchange and correlation. Additional IXS results obtained for K and Rb are addressed briefly. 



PACS numbers: 71.45.Gm, 62.50.-p, 78.70.Ck, 71.20.-b 

Sodium, at ambient conditions, is one of the best manifes- 
tations of a "simple" or nearly-free-electron (NFE) metal [l]. 
It is characterized by a single s-type valence electron, weak 
interaction between the conduction electrons and the atomic 
cores (electron-ion interaction), and conduction band states 
of orbital character Na crystallizes in the high-symmetry 
body-centered cubic (bcc) structure at pressures up to 65 GPa, 
where it transforms to face-centered cubic (fee) [2J. The 
properties of Na change fundamentally under pressure in the 
megabar pressure range, where a series of phase transitions 
into lower-symmetry crystal structures has been predicted 13J 
and observed |4-6|, accompanied by marked changes in its 
optical properties |4J |6l |2l and culminating in the formation 
of a non-metallic, visually transparent phase at ~200 GPa |8|. 
A central question is how the transformation from a simple 
metal to a semiconductor progresses, not only in Na, but also 
in other metals such as Li, which was reported to become 
semiconducting above 70 GPa |9|. As for Na, does it remain 
NFE-like in its bcc and fee phases up to ~100 GPa |7 1 so that 
the non-NFE behavior starts only with the transitions into the 
lower-symmetry phases above 105 GPa, or are there signifi- 
cant precursors at lower pressure? 

To provide an answer, we measured and calculated the pres- 
sure dependence of Na plasmon energies. Plasmon excita- 
tions provide information on the collective electronic excita- 
tions in the form of longitudinal charge density waves at fi- 
nite wavevector, and they determine the optical response of a 
metal, specifically the plasma reflection edge. Plasmons have 
been studied for many years by electron energy loss spec- 
troscopy (EELS) at zero pressure (see for example ifTOl fTTl 
and references therein), but this technique is not suitable for 
samples enclosed in high-pressure cells. Mao et al. fY2\ 
have demonstrated the possibility of measuring plasmon ex- 
citations in Na under pressure using inelastic x-ray scattering 
(IXS), and they found their experimental results up to 2.7 GPa 
to be in agreement with theoretical predictions. 

We report here detailed IXS results on the plasmon energy 
dispersion in Na under pressure up to 43 GPa, corresponding 
to a 2.6-fold increase in density. Our results evidence a sig- 
nificant departure from the predictions for a NFE metal. In 
order to explain this discrepancy between theory and experi- 



ment, we also present plasmon energy calculations based on 
first-principles electronic band structures and a quasi-classical 
plasmon model after Paasch and Grigoryan (PG model) 1 1 3 1 . 
These calculations reconcile experiment and theory and allow 
us to identify changes in the electron-ion interaction as the 
dominant effect, rather than changes in the electron-electron 
interactions. Some experimental results are also reported for 
K and Rb. 

The theoretical description of plasmons in simple metals is 
well established. The most commonly used approach starts 
from the free-electron (EE) gas and uses the Random Phase 
Approximation (RPA) |14|. The plasmon energy dispersion 
Ep{q) is then given by Ep{q) - tiUp + ^ aq with the plasma 
frequency a)p - -y/ne^/eoejm and the dispersion coefficient 
ffPE = ^Ef/Ticijp, where q is the plasmon momentum, m the 
electron mass, n the electron density, a dielectric constant 
describing the polarizability of the ionic cores (e, = 1 for the 
free-electron gas), and Ef - (h^ /2m)(37T^n)^^^ is the Fermi 
energy. This relatively simple model works reasonably well 
for simple metals such as Na and Al at ambient conditions 
ifTOirm . but experimental dispersion coefficients a tend to be 
lower than the theoretical values, which has been attributed 
to electron exchange and correlation effects |11|. Both the 
plasma frequency ojp and the dispersion coefficient a depend 
on the electron density and can thus be tuned by the applica- 
tion of pressure; both are expected to increase with increasing 
pressure. The stability of bcc sodium over a large pressure 
range of 0-65 GPa permits to generate an up to threefold in- 
crease in (electron) density without a structural transition. 

IXS experiments were performed on beamline ID 16 at the 
ESRF, Grenoble. Silicon crystal monochromators were used 
to monochromatize the incident beam and to analyze the scat- 
tered radiation. The incident x-ray beam with a photon energy 
of 9.877 keV was focussed onto the polycrystalline sample 
in a high-pressure cell with a spot diameter of 100-200 fim, 
depending on the sample size. IXS spectra of the samples 
at room temperature were recorded in energy-scanning mode 
with an overall spectral resolution of 0.6 eV and a momentum 
resolution of 0.4 nm"'. In most of the experiments, the sam- 
ples were pressurized in diamond anvil cells (DACs). Rhe- 
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FIG. 1. IXS spectra of polycrystalline sodium pressurized in a dia- 
mond anvil cell, (a) Energy transfer spectra for a momentum transfer 
of ^ = 5 nm"' and pressures of 3.5-43 GPa. (b) Energy transfer spec- 
tra of Na at 16 GPa and momentum transfers of 2-9 nm"'. Vertical 
offsets are added for clarity. 



nium and stainless steel gaskets were used with initial thick- 
nesses of 50-100 jum and hole diameters of 150-200 //m. The 
loading of distilled Na, K, and Rb into the pressure cells was 
carried out in an argon atmosphere. Because of the softness 
of these metals no pressure transmitting medium was added. 
Pressures were determined with the iTiby method |15| or by 
measuring an x-ray powder diffraction of the sample and us- 
ing its known equation of state fTl. In the experiments us- 
ing a DAC, the incoming x-ray beam passed through one dia- 
mond anvil onto the sample, and the scattered radiation was 
collected through the opposing anvil (thickness ~1.5 mm). 
Despite the small sample scattering volume and the relatively 
low diamond transmittance of 30% for 10-keV x-rays, high- 
quality spectra could be collected in ~2 hours. The experi- 
ments on Na at 1 and 7 GPa were performed using a high- 
pressure cell equipped with sintered-diamond anvils and a 
beryllium gasket. Here, the x-rays passed through the Be gas- 
ket. 

Figure [T] shows IXS spectra of polycrystalline Na under 
pressure. A single excitation peak is observed, which is at- 
tributed to the bulk plasmon in Na, based on the predictions 
for a NFE metal with the electron density of Na and also 
by comparison with previous EELS results fTOl. The rising 
background in the spectra of Fig. [T] in particular at large q, is 
due to plasmon and interband excitations in the diamond anvil 
through which the scattered radiation is detected |16|. 

Plasmon energies and linewidths were determined by fitting 
the spectra with a Gaussian peak for the plasmon line and a 
polynomial background. As all experiments were performed 
on polycrystalline samples, the reported plasmon energies are 
directional averages. The Na plasmon energies increase with 
increasing pressure as illustrated in Fig. |2] for <? = 5 nm"'. 
This is in qualitative agreement with the NFE picture, where 
the plasmon energies scale with the electron density. 

Figure|2]also shows results for K and Rb. These two metals 
clearly do not follow the expectations for NFE metals. More- 



FIG. 2. Plasmon energies in Na, K, and Rb at iji = 5 nm"' as a func- 
tion of pressure. Ambient-pressure EELS results B lOl are indicated 
by stars. For the Na data, different symbols distinguish different sam- 
ple loadings and pressure cells. Lines are guides to the eye. 



over, their plasmon linewidths increased and thek plasmon 
intensities decreased rapidly with increasing pressure. These 
effects are attributed to the pressure-driven s-d hybridization 
of conduction band states, as is also evident from the optical 
response of K and Rb under pressure ifTTlfTSl (see also [19|). 
We will therefore focus on Na that could be studied over the 
largest pressure range and which, at ambient conditions, is one 
of the best manifestations of a nearly-free-electron metal. 

Figure |3ja) shows the measured plasmon dispersion rela- 
tions of Na at several pressures. The plasmon energies are 
plotted versus (f' because of the anticipated parabolic dis- 
persion relation, see the Ep{q) relation given above. At the 
lowest pressure, 1 GPa, the Na plasmon dispersion measured 
here is indeed very close to parabolic. Results of an ambient- 
pressure EELS study |10| are included in Fig.[3|a) for com- 
parison. In the low-g' region, the IXS and EELS data are rea- 
sonably consistent, but the EELS results exhibit some devia- 
tion from a parabolic dispersion. The offset, at low q, between 
the IXS and EELS data is largely due to the pressure applied 
in the IXS experiment. The 1-GPa plasmon dispersion mea- 
sured by IXS is described best by £(0) = hiOp - 5.82(2) eV 
and a - 0.26(1). As noted before [ 101 , the experimental 
dispersion parameter a is lower than the FE/RPA value of 
QJRPA = 0.35. Towards higher pressures, the measured plas- 
mon dispersion remains approximately parabolic, deviations 
being most notable at 16 GPa. 

Figure [3|b) shows that also the plasmon linewidth is 
strongly pressure dependent. This effect is tentatively at- 
tributed to a reduction of the plasmon Ufetime due to decays 
involving electron-hole excitations. A detailed study by, e.g., 
time-dependent density functional theory could be a subject 
for future studies. For the remainder of this paper we focus on 
the plasmon energies, assuming that self-energy effects on the 
plasmon frequency can be neglected. 

The IXS results on the plasmon energies of Na as a func- 
tion of momentum and pressure allow us to test the validity 
of the FE/RPA description of Na and to assess the relative 
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FIG. 3. (a) Plasmon dispersions E{(f-) of polycrystalline Na as a 
function of pressure. Ambient-pressure electron-energy-loss spec- 
troscopy results by vom Felde etal. 1 10| are indicated by small open 
symbols, (b) Plasmon linewidth of Na versus pressure. 



FIG. 4. Experimental pressure dependence of the directionally- 
averaged plasmon energy in Na at ^ = 5 nm"' (large symbols) and 
results of the free-electron (FE) gas model, EE model with core polar- 
ization, and the PG model (solid line). The star indicates the ambient- 
pressure EELS result 1101 . 



importance of band structure effects, core polarizability, and 
exchange-correlation. Using Ep(q) as given above and the ex- 
perimental equation of state of Na [T|, the pressure depen- 
dence of the q = 5 nm ' plasmon was calculated within the 
FE/RPA framework as shown in Fig. |4] This FE/RPA esti- 
mate is significantly higher in energy than the experimental 
values, and the deviation increases with increasing pressure. 
Inclusion of the core polarization (fj - 1.16, Ref. Il20l ) leads 
to a good agreement with the experiment at low pressure, but 
a major deviation between theory and experiment remains at 
high pressure. 

The important observation here is the striking increase in 
the deviation between theory and experiment with increas- 
ing pressure. In previous work on other metals, deviations 
from the FE/RPA predictions were discussed in relation to ex- 
change and correlation effects lfT0ll2TII . and a number of exten- 
sions of the RPA were proposed in this spirit (see IfTOl 12111221 
and references therein). Their main effect is to reduce the 
plasmon dispersion coefficient a, and their inclusion can im- 
prove the agreement between theory and experiment. How- 
ever, exchange-correlation effects decrease with increasing 
electron density, and they can thus be excluded as the origin 
of the deviation between the 'FE/RPA + core polarization' re- 
sults and the experimental data in Fig. [4] More recent theoreti- 
cal studies have emphasized the importance of band-structure 
effects on the plasmon properties, regarding both the energy 
dispersion |13| and the plasmon linewidth 119]. The present 
IXS results offer an opportunity to test these proposals. 

The PG plasmon model lfT3l adopted here is an extension 
of the FE/RPA approach. This classical model is not ex- 
pected to describe the plasmon propeities as accurately as, 
e.g., time-dependent density functional theory (DFT) lfT9l . 
but it allows us better to understand the underlying physics. 
The electronic structure of the metal is described here by a 
single isotropic conduction band with a quartic dispersion, 
E(k) = E2k^ + E^k^ . The E^k^ term accounts for deviations 
from the parabolic band shape of the free-electron gas. The 



plasma frequency Up and the plasmon dispersion coefficient 
a can then be determined from Eo, E4 and the Fermi energy 
Ef as described in detail in 1 13 1. 

Electronic structure calculations of Na were performed in 
the framework of first-principles DFT, using the full-potential 
L/APW-Hlo method 1 23 , 24 1 and the Generalized Gradient Ap- 
proximation II25I as implemented in the WIEN2K code ll26l . 
The electronic band structure of bcc Na was calculated for a 
series of volumes corresponding to the pressure range of 0- 
50 GPa Ii27 |. The coefficients E2 and E4 were determined for 
three directions in the Brillouin zone, i.e. along [^00], [^^0], 
and [^^^], and then averaged |28 1. The only adjustable param- 
eter in the calculation of the plasma frequency and the disper- 
sion coefficient is the static dielectric constant that accounts 
for the polarizability of the ionic cores. A value of - 1.16 
as determined for Na in a DFT-based study Il20l was adopted 
here, and a possible density dependence was neglected. 

Figure|4]shows the results of the PG model with input from 
our band structure calculations for the q - 5 nm ' plasmon 
in Na [281 • The computed results are in excellent agreement 
with the IXS data. To trace the source of the difference be- 
tween the FE/RPA results and those of the PG model, Table |l] 
summarizes the plasma frequency Hojp and the dispersion co- 



TABLE I. Experimental plasma frequencies cjp and dispersion coef- 
ficients a of Na as a function of pressure and comparison with theo- 
retical results. The experimental values were determined from linear 
fits to the E(q~) data as shown in Eig.jsja). 
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efficient a for the two models at selected pressures. At 1 GPa, 
the results of the NFE/RPA and the PG model are very similar, 
and the calculated plasma frequency agrees well with the ex- 
periment. As noted before |10|, the calculated values of a are 
20-30% larger than in the experiment, and this is probably due 
to exchange-correlation effects not included here. The effect 
of the non-parabolic contribution (£4 < 0) is to reduce both 
tiojp and a. In Na, this effect is very small near ambient pres- 
sure, confirming the analysis by Paasch and Grigoryan |fT3l . 
With increasing pressure, however, |£'4| increases and leads to 
substantial corrections |28|. At 40 GPa, it reduces the plasma 
frequency by 8% and the dispersion coefficient by 35% com- 
pared to NFE/RPA. Table |l] also shows that a decreases with 
increasing pressure in the PG model, in contrast to the free- 
electron behavior. We would like to emphasize that it is the 
non-parabolic contribution that causes the renormalization of 
the plasmon energies, even though the corrections to the band 
energies are less than 3% of the band width |28 1. 

As for the physical origin of these corrections, the band 
structure calculations show that pressure causes the band gaps 
at the N, P, and H points of the Brillouin zone to grow rela- 
tive to the width of the conduction band ESI . This evidences 
a strengthening of the interaction between the valence elec- 
trons and the ionic cores (electron-ion interaction). In other 
words, Na becomes increasingly less free-electron like un- 
der compression. The quartic correction (£4 < 0) corre- 
sponds to a lowering of the conduction band energies near the 
Brillouin zone boundary in comparison to the free-electron 
case, and this distortion of the band structure |28 | leads to 
a reduction of the Fermi velocity, Vf - {\IK)dE{k)ldk\kp - 
(2kf/K)(E2 + lEi^k^p), which is the key physical quantity that 
determines the plasmon dispersion. 

In summary, we have performed inelastic x-ray scattering 
experiments to determine the effect of pressure on the plasmon 
excitations in the 'simple metal' sodium. While Na is consid- 
ered one of the best manifestations of a nearly-free-electron 
metal at ambient conditions, our results evidence substantial 
and increasing deviations from the behavior of a NFE metal at 
high pressure up to 43 GPa. This can be seen as an early pre- 
cursor of the fundamental changes in the electronic structure 
of Na at megabar pressures. The deviation from NFE behavior 
can in part be attributed to the polarizability of the ionic cores, 
but to a larger extent it is caused by pressure-induced changes 
in the electronic band structure. Plasmon energies determined 
on the basis of electronic band structure calculations and the 
quasi-classical PG model are in excellent agreement with the 
IXS results. They show that the electron-ion interaction in 
Na increases with pressure and leads to the renormalization 
of the plasmon energies via a modification of the electronic 
band structure. As for the heavy alkali metals under pressure, 
band structure effects can be expected to be even more impor- 
tant due to the pressure-driven hybridization of the valence s 
orbitals with d states, as discussed before |13|. We observed 
a weak pressure dependence of the plasmon frequencies of 
K and Rb combined with fast broadening of their plasmon 
resonances under pressure. These results may aid the inter- 



pretation of plasmon dispersions of the heavy alkali metals at 
ambient pressure. As for bulk Li metal, the quite detailed pre- 
dictions on the collective electronic response under pressure 
11291 |30l still await a related experimental investigation. 
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